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Abstract — Optimum Bayes estimator for General Gaussian 
Distributed (GGD) data in wavelet is provided. The GGD 
distribution describes a wide class of signals including natural 
images. A wavelet thresholding method for image denoising is 
proposed. Interestingly, we show that the Bayes estimator for this 
class of signals is well estimated by a thresholding approach. This 
result analytically confirms the importance of thresholding for 
noisy GGD signals. We provide the optimum soft thresholding 
value that mimics the behavior of the Bayes estimator and 
minimizes the resulting error. 

The value of the threshold in BayesShrink, which is one of 
the most used and efficient soft thresholding methods, has been 
provided heuristically in the literature. Our proposed method, 
denoted by Rigorous BayesShrink (R-BayesShrink), explains 
the theory of BayesShrink threshold and proves its optimality 
for a subclass of GDD signals. R-BayesShrink improves and 
generalizes the existing BayesShrink for the class of GGD signals. 
VV^hile the BayesShrink threshold is independent from the wavelet 
coefficient distribution and is just a function of noise and noiseless 
signal variance, our method adapts to the distribution of wavelet 
coefficients of each scale. It is shown that BayesShrink is a 
special case of our method when shape parameter in GGD is 
one or signal follows Laplace distribution. Our simulation results 
confirm the optimality of R-BayesShrink in GGD denoising with 
regards to Peak Signal to Noise Ratio (PSNR) and Structural 
Similarity (SSIM) index. 

Index Terms — Bayesian estimation, wavelet shrinkage and 
denoising, soft thresholding 

I. Introduction 

Data denoising via wavelet thresholding is a simple yet 
powerful method that has been studied and used in various 
research areas ranging from communications to biomedical 
signal analysis. Wavelet shrinkage is based on rejecting those 
wavelet coefficients that are smaller than a certain value 
and keeping the remaining coefficients. Thus, the problem of 
removing noise from a set of observed data is transformed 
into finding a proper threshold for the data coefficients. The 
wavelet denoising method consist of three steps: 1) applying 
the wavelet transform to data and calculating the wavelet 
coefficients; 2) applying hard or soft thresholding method; and 
3) calculating the denoised signal by using the inverse wavelet 
transform. The pioneer shrinkage methods, such as VisuShrink 
and SureShrink, propose thresholds that are functions of the 
noise variance and the data length Q-El- Over the past fifteen 
years, several thresholding approaches such as ||4l-||7l have 
been developed. These methods provide optimum thresholds 
by focusing on certain properties of the noise-free signal, and 
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they are proposed for particular applications of interest. In 
wavelet shrinkage methods, the noise component is usually 
assumed to have a Gaussian distribution, whereas various 
distributions have been considered to model the noiseless 
coefficients. Some of the signal classes are as follows: family 
of Pearsons distributions ID, Laplace distribution [9|, family 
of SaS distributions fTOl, or Gauss- Markov random field 
model ifTTIl . It is important to note that, in general, the 
distribution of wavelet coefficients changes from scale to scale. 
Two approaches have been used by different groups to deal 
with this change. One approach is to use a fixed model in all 
scales which decreases the modeling and denoising accuracy. 
The second solution is to fit the considered distribution of the 
wavelet coefficients at each scale by optimizing the parameter 
in that scale 0, ifTOl . ifTTIl . However, it is known that wavelet 
transform is not a complete decorrelator. Secondary properties 
in wavelet transform such as clustering and persistence result 
in Large/small values of wavelet coefficients to propagates 
across the different scales. To model these properties, |[T3ll 
and [14J have jointly used Gaussian model, which can eas- 
ily capture the correlation between the wavelet coefficients. 
Nevertheless, this model is unrealistic since the wavelet co- 
efficients have heavy tailed distribution and in addition, the 
proposed models are unable to include frequent sparsity and 
compression properties of wavelet coefficients. Consequently, 
wavelet based Hidden Markov models (HMM) that capture 
the dependencies between the coefficients with non-Gaussian 
distribution have been proposed |[T2l . Although these methods 
outperform the thresholding methods, they are computationally 
expensive and time consuming which make them inefficient 
when speed and time are of great importance. In these cases 
we compromise between complexity and performance and lean 
towards thresholding algorithms. 

In this paper we aim to find the Bayesian estimate of 
the noiseless signal from its observed noisy wavelet coeffi- 
cients for General Gaussian Distributed (GGD) signals. We 
assume that the subbands of wavelet transformation of the 
images follow General Gaussian distribution, in which the 
shape parameter is often between 0.5 and 1. This distribution 
is highly accepted in the image processing community for 
wavelet coefficients of image subbands llT5ll - |[T9l . Our results 
show that the optimum Bayes estimator in this class of signals 
behaves similarly to a soft thresholding method. Therefore, 
it confirms the importance and efficiency of the existing soft 
thresholding approaches for image denoising purposes. One of 
the well known and most efficient image thresholding methods 
is BayesShrink ||4|, which provides the threshold based on 
minimizing the Bayes risk heuristically. The threshold value 
proposed in BayesShrink is a function of noise variance and 
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variance of noiseless data and is independent from the shape 
parameter. However, since distribution of wavelet coefficients 
vary in each scale, using a threshold value which is not 
adapted for each scale is not very accurate. Here we propose a 
generalized threshold value which adapts to distribution shape 
of coefficients in each scale. Simulation results confirm that the 
newly proposed method denoted by Rigorous BayesShrink(R- 
BayesShrink) outperforms the BayesShrink itself. 

The paper is organized as follows. In Section [HI the de- 
noising problem is formulated. In Section Jill we calculate the 
Bayes estimator for GGD signals. The proposed thresholding 
algorithm based on the Bayes estimator is provided in Section 
HVl Section |Vl contains the simulation results, and finally in 
Section |Vl] conclusions are drawn. 

II. Problem Formulation 

We denote the noise-free data of length N with vector = 
[^1, • • • jVn]^- The noise-free data is corrupted by an additive 
white Gaussian random process = [wi^ - ■ ■ jWn]^ with 
zero mean and variance of a^. Therefore, the observed data 
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which can be expressed in terms of a desired orthonormal 
wavelet basis in the following form 



(2) 



where Oi is the ith noisy wavelet coefficients, Oi is the ith 
noiseless coefficient and ViS are the noise wavelet coefficients. 
Since orthonormal bases are used, ViS are also additive white 
Gaussian random variables with the same mean and variance 
as the noise samples (wi). 

We assume that the noiseless coefficients (0) have zero mean 
General Gaussian Distribution EI, 117611 - |[T9ll : 



f{0) = C{ay,l3)exp{-[a{ay, 



(3) 



where f3 is the shape parameter, a{ay^f3) = ay [^ry] , 
C{cTy^ P) = and r{x) = t-^e-'dt. 



When P is one, GGD is the same as Laplacian distribution 
and when /3 is two, GGD describes the Gaussian distribution. 
The objective of this work is to provide the optimum soft 
threshold, T, that minimizes the Bayes risk. This question 
is addressed by BayesShrink approach heuristically and its 
threshold is provided by numerical analysis. Our goal is to 
search for the optimum threshold analytically. 

As a reminder, in both hard Q and soft © thresholding, 
small absolute values of are mapped to zero |31. The hard 
thresholding keeps the values larger than the threshold. So the 
map is along = after that threshold. 



if 0{i) < T, 
9{i) if otherwise. 



(4) 



While the soft thresholding map is off this diagonal line by 
the amount of the threshold. 



eT{i) = sign{0{i)mi)\-T)+ (5) 
where (x)+ = max(x,0). 

III. Bayesian Estimate of GGD Signals 

Here we calculate the general Bayes estimate structure for 
the GGD signals. The Bayes estimate of the desired noise- 
free parameter 0, denoted by 0, minimizes the following mean 
square error (MSE) 

Ob = d^Tg^min E{0 -Of (6) 

This estimator is the equivalent of the mean of the posterior 
distribution that is an unbiased least- squares estimate of 0, 
given measurement f?! 

Due to the Gaussian distribution of the additive noise, we have 



fe\em = fv{0-0) 
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and thus Q can be written as 

^ _ j0M0-0)fsi0)dO 

jM0-0)f5{O)d0 

For (3 = 2 (the Gaussian case) this Bayes estimator will 
become 

(9) 



where cr| is the variances of noiseless data. 

While equation ^ cannot be solved for all other values of 
/3, it will have the following analytical form for P = 1 (the 
Laplacian case) 

V2e 
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in which Q{x) = e 2 dt. Figure [T] shows this 

Bayes estimator, Ob, for different values of ciy while cfw is 
constant and 

IV. Mean Square Error(MSE) Soft Threshold 
Fitting for the Bayes Estimator 

We studied the behavior of the Bayes estimator of noisy 
GGD signals, in dl]), as the noise variance, the noise-free 
variance and the shape parameter vary. Our analysis confirms 
that the estimator behaves similarly to a soft thresholding 
algorithm [5] as will be described. Figures [T] and [2l show 
examples of variation of this estimator. 

Note that a soft thresholding map is an odd function of 
the noisy coefficient and is zero up to a particular point T 
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Fig. 2. Bayes estimator curves for different noise variance cr xu while noiseless 
data variance Gy is five and shape parameter /3 is one. 



and then follows a line with slope of one. Our observation, 
as depicted in these two figures, confirms that the Bayes 
estimate is an odd function of the noisy coefficients that is 
zero at zero, grows nonlinearly with a small amount to a 
particular point and after that it asymptotically follows a line 
for larger wavelet coefficients. In Figure [3] we show the values 
of this slope as the shape parameter /3 varies. As the figure 
shows, the slope is extremely close to one for when the same 
parameter is between zero and one. This behavior mimics the 
soft thresholding. It is known that in natural images we deal 
only with the shape parameters between 0.5 to 1 llH, |[T6l - 
lfT9l . Therefore, a soft thresholding estimate of the Bayes 
estimate can very well mimic the behavior of this estimator 0. 
To find the optimum thresholding mapping, we minimize the 
following least square error between the Bayes estimate in ([8]) 

^Note that for the /5's out of this range, we still can use weighted version 
of the coefficients based on Figure [3] and then provide a thresholding. 



and the soft thresholding by value T in ([5]): 

/ + 00 
!{^){BB-QTfAB (11) 
-CO 

Figure H] shows these threshold values found for a range of 
noise variance, signal variance, and shape parameter. Our 
study, as depicted in this figure, resulted in the following 
choice of a closed form representation that best describes this 
threshold value as a function of noise variance (a^y), noiseless 
variance {oy), and shape parameter 

where parameters a, h\, 62 and bs{beta) can be calculated 
through another least square optimization. The following op- 
timum values have been found by solving this optimization 
problem over a feasible range of noise variance, signal vari- 
ance, and shape parameter: 

a = 1, 61 = 62 = 1 + V^, bs = (13) 

The above parameters lead to the following optimum soft 
threshold that best mimics the Bayes estimator 

T^ = ^a.(^)^ (14) 
VP 

We denote this threshold as Rigorous BayesShrink (R- 
BayesShrink). 

A. BayesShrink and Rigorous BayesShrink 

We showed in the previous sections that the Bayes estimator 
for GGD signals behaves very much like a soft thresholding 
algorithm. Therefore, we have shown rigorously that soft 
thresholding can indeed capture the behavior of Bayes es- 
timator in this scenario. BayesShrink is one of the most 
efficient existing methods of wavelet thresholding for the 
purpose of image denoising lH. The method searches for the 
optimum soft threshold that minimizes the MSE and maps the 
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observation to a class of thresholds in the form of cr^ . Its 
numerical analysis finds the following threshold 



(15) 



Nevertheless, if we set the value of P to one in ([T4l) . we 
get the same threshold of the BayesShrink. This confirms 
and rigorously proves the importance of BayesShrink for the 
Laplacian case. However, as the BayesShrink threshold ignores 
the role of P of the considered GGD model, it is expected 
that R-BaysShrink outperform the BayesShrink. Examples 
in our simulation results confirm that R-B ay Shrink indeed 
outperforms BayesShrink. 

B. Least Square Estimator(LSE) vs. Rigorous BayesShrink 

The optimization for soft threshold mapping in R- 
BayesShrink is based on minimizing the MSE in (fTTI) . In 
an earlier work, we have optimized this mapping by LSE 
minimization II2TII 

/+00 
\eB-0T?de (16) 
-00 

It was shown that using the structure in ([T2I) . the optimum 
threshold is 



Tlseb = 



cr' 



(17) 



We denote this threshold as LSE Based (LSEB). 

On the other hand, the Maximum a posterior (MAP) esti- 
mator for the Laplace distribution (/3 = 1) is 

Omap = argm_ax[log(/,(^ - 0)) + log(/,-(^))] (18) 
e 

and since the additive noise follows a Gaussian distribution, 
the MAP estimate is 

{0 - Of 



Omap argmax[- 



log{awV2) 



(19) 



Solving this optimization, leads to the following soft threshold 
value 



Tmap 



ay 



(20) 



This value is the same as the LSE based threshold mapping 
in ([TtI) when the P is set to one. This confirms that for the 
Laplacian distibution, the LSE based estimator is the same as 
the MAP estimator! Note that this threshold value is the same 
as BayesShrink (R-BayesShrink for /3 = 1) with an extra term 
V2. 

C. Parameter Estimation of GGD 

The proposed R-BayesShrink threshold ([T4l) has three pa- 
rameters: 

1) noise variance {(Jw)'- Median Absolute Deviation (MAD) 
estimate is the standard noise variance estimator d, 1201 : 

Median{\e\) 



O'MAD 



0.6745 



-,6> G subband HHi (21) 



TABLE I 

Shape parameter ((3) estimated from different noisy data with 
different snr's using mme. 



/5 


0.4 


0.6 


0.8 


1.0 


SNR=15 


0.45 


0.69 


0.83 


1.01 


SNR=30 


0.41 


0.68 


0.83 


1.03 



where HHi is the high frequency subband. 

2) Noiseless data variance (cry): Since the data is indepen- 
dent from noise, the variance of the noisy data is 



2 2,2 
ay = ay^a^ 



On the other hand, variance of the noisy data can be estimated 
by: 



(22) 



yieN{k) 



where M is the number of coefficients in scale N(k) and rriy 
is the mean in the same scale my = Using the noise 

variance estimate in ([2T]) and noisy data variance estimate in 
the estimate of noiseless data variance is 



-I 



(23) 



3) Shape parameter (/3): There are three methods available 
for P estimation in Generalized Gaussian Distributed signals: 
A) Moment Matching Estimator (MME), which uses moments 
of GGD signals; B) Entropy matching Estimator (EME) that 
relies on matching the entropy of the GGD modeled distri- 
bution with that of a set of empirical data; and C) Maximum 
Likelihood Estimator (MLE), which finds the shape parameter 
which maximizes the likelihood ll22l . 

MME is an accurate and still a simple method that is useful 
when p is between 0.18 and 2 which covers the range for 
images [|22l. Therefore, we use MME for GGD parameter 
estimation. 

To estimate /3, MME uses the second and fourth moments 
of the subbands. For the Kurtosis of GGD, i^y\y-eN{k)^ we 
have (H 

1/^22 Q 4 , / 2 2 \2^'^^)^'^f \ 

(24) 

and i^y\y-eN{k) itself can be estimated from the observed noisy 
data {yi\yi G N{k),i == 1,...,M} as follows ll23ll : 



f^y\yieN{k) 



(25) 



M 



(M - l)al 

Where rriy is the mean value of available i/iS (rriy 
M is the number of coefficients in k^^ scale and ay is from 
(ESI)). 

Having the estimate of ny\y.^N{k) from (l25l) , estimate of 
noise variance from (I2TI) . and ay estimate from (1231) , we can 
estimate {3 from (l24l) . Table [J shows some of the estimated 
/3's from the given noisy data. It is evident that the estimated 
values are reasonably accurate. 
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V. Simulation Results 

Here we compare the performance of our proposed algo- 
rithm and BayesShrink for three sets of data: 

1) Synthetic data, which is a set of GGD data with a known 
shape parameter. 

2) Natural images. 

3) Computed Tomographic (CT) images. 

In all the tests, we use Haar wavelets in 5 levels. Each test is 
run one hundred times and image sizes are 512 x 512. 

To measure the quality of the denoised images. Peak Signal 
to Noise Ratio (PSNR) is used which is defined as follows: 

(2^ - 1)2 

PSNR=lQlog,o ^ ^g^^ (26) 

where B is the number of bits in the images and MSE is the 
Mean Square Error (MSE). For natural images B = 8 and for 
CT phantom images B = 16. 

Table [III shows the PSNR of denoised synthetic data with 
BayesShrink in ([15]), LSEB in ([17]), and R-BayesShrink in 
([T4l) . Based on this table, R-BayesShrink outperforms the other 
two methods and as was expected, for /3 = 1, R-BayesShrink 
and BayesShrink performance are the same. 
Since R-BayesShrink performs better than LSEB, we elimi- 
nate the comparison of LSEB in the following examples. 

For the natural images (CameraMan, Lena, Mandrill, pep- 
pers, goldhill, and boatshown) the denoising results are given 
in Table Uni 

Table Hill shows that R-BayesShrink has better PSNR 
values which shows its better performance in comparison 
with BayesShrink. Figure [5] shows the denoised images 
with BayesShrink and R-BayesShrink. It is evident that R- 
BayesShrink preserves the edges better than BayesShrink. To 
quantitatively show this property. Structural Similarity (SSIM) 
index [24] is used, and the results are given in TabelHV] SSIM 
index compares local patterns of pixel intensities which are 
normalized for luminance and contrast. 

The third test is done with CT images with different x- 
ray doses. In these images the x-ray tube voltage is kept at 
120kV, and its current (mA) is changed. The quality of the 
CT images is a function of x-ray tube current and voltage. 
Decreasing the tube current causes a decrease in the number of 



TABLE II 

PSNR OF SYNTHETIC DATA DENOISED WITH BAYESSHRINK, LSEB, AND 

R-BayesShrink. 



SNR 1 5 1 10 


15 1 20 1 25 1 30 


/3= 


0.1 


Bayesshrink 


35.22 


37.29 


45.36 


49.77 


52.22 


59.19 


LSEB 


36.63 


38.24 


45.82 


50.02 


52.35 


59.27 


R-B ayes Shrink 


37.77 


44.97 


49.23 


53.10 


55.17 


63.10 




0.3 


Bayesshrink 


33.15 


37.56 


41.68 


44.63 


52.46 


56.85 


LSEB 


33.83 


37.99 


41.91 


44.75 


52.51 


56.87 


R-B ayes Shrink 


34.08 


38.47 


42.41 


45.18 


52.77 


57.04 


/3= 


0.5 


Bayesshrink 


27.74 


30.73 


38.45 


40.43 


43.92 


47.40 


LSEB 


27.97 


30.92 


38.56 


40.47 


43.95 


47.40 


R-B ayes Shrink 


27.97 


30.98 


38.65 


40.50 


43.97 


47.41 


13= 


0.8 


Bayesshrink 


21.03 


24.52 


33.28 


33.84 


39.28 


45.07 


LSEB 


20.98 


24.55 


33.31 


33.85 


39.29 


45.08 


R-B ayes Shrink 


21.05 


24.56 


33.31 


33.85 


39.29 


45.10 


/?=i 


Bayesshrink 


18.99 


23.26 


29.36 


33.71 


40.18 


44.58 


LSEB 


18.80 


23.24 


29.35 


33.70 


40.17 


44.56 


R-B ayes Shrink 


18.99 


23.26 


29.36 


33.71 


40.18 


44.58 



TABLE III 

PSNR OF DENOISED IMAGES WITH BAYESSHRINK AND THE PROPOSED 
THRESHOLD (R-BAYESSHRINK). 



SNR 


5 


10 


15 


20 


25 


30 


CameraMan 


BayesShrink 


20.99 


23.14 


25.92 


29.25 


32.99 


37.18 


R-B ayes Shrink 


21.08 


23.24 


26.09 


29.41 


33.20 


37.28 


Lena 


BayesShrink 


23.89 


25.89 


28.14 


30.78 


33.82 


37.35 


R-B ayes Shrink 


23.94 


25.96 


28.27 


30.91 


34.00 


37.54 


Mandrill 


BayesShrink 


20.04 


21.52 


23.94 


27.40 


31.46 


35.98 


R-B ayes Shrink 


20.07 


21.57 


23.94 


27.41 


31.48 


36.00 


GoldHill 


BayesShrink 


24.00 


25.68 


27.67 


30.38 


33.37 


37.42 


R-B ayes Shrink 


24.04 


25.71 


27.73 


30.38 


33.54 


37.48 


Peppers 


BayesShrink 


20.78 


23.29 


26.34 


29.44 


32.94 


36.87 


R-B ayes Shrink 


20.83 


23.37 


26.41 


29.76 


33.16 


37.00 


Boat 


BayesShrink 


22.48 


24.35 


26.68 


29.42 


32.64 


36.45 


R-B ayes Shrink 


22.52 


24.43 


26.77 


29.55 


32.80 


36.57 



Fig. 5. Denoised images with BayesShrink and R-BayesShrink. Top: Noisy images with SNR=15, middle: denoised with BayesShrink and bottom: denoised 
by R-BayesShrink. 



TABLE IV 

SSIM OF DENOISED IMAGES WITH BAYESSHRINK AND THE PROPOSED 
THRESHOLD (R-BAYESSHRINK). 



SNR 


5 


10 


15 


20 


25 


30 


CameraMan 


BayesShrink 


0.13 


0.30 


0.54 


0.77 


0.90 


0.97 


R-BayesShrink 


0.14 


0.33 


0.59 


0.83 


0.96 


0.99 


Lena 


BayesShrink 


0.33 


0.60 


0.81 


0.94 


0.98 


0.99 


R-BayesShrink 


0.35 


0.60 


0.83 


0.95 


0.99 


0.99 


Mandrill 


BayesShrink 


0.35 


0.74 


0.93 


0.993 


0.998 


0.999 


R-BayesShrink 


0.45 


0.79 


0.95 


0.996 


0.999 


0.999 


GoldHill 


BayesShrink 


0.32 


0.63 


0.85 


0.95 


0.993 


0.999 


R-BayesShrink 


0.35 


0.65 


0.89 


0.97 


0.996 


0.999 


Peppers 


BayesShrink 


0.14 


0.33 


0.58 


0.79 


0.92 


0.97 


R-BayesShrink 


0.16 


0.35 


0.63 


0.85 


0.96 


0.99 


Boat 


BayesShrink 


0.34 


0.64 


0.85 


0.95 


0.990 


0.998 


R-BayesShrink 


0.37 


0.65 


0.87 


0.98 


0.995 


0.999 



photons hitting the detectors and thereby decreases the quaUty 
and SNR. Here we test the images with 8 different currents 
from 10mA to 500mA. To calculate the PSNR we use the 
images captured by 500mA as the reference and compare the 



TABLE V 

Noise standard deviation after and before denoising with 
BayesShrink AND R-BayesShrink which is calculated based on 

SIGNAL VARIATION IN A SMOOTH REGION (THIS IS HOW RADIOLOGISTS 
CALCULATE THE AMOUNT OF NOISE). 



Tube Cur. (mA) 


10 


20 


50 


100 


200 


300 


400 


500 


Original 


79.70 


59.50 


29.10 


18.80 


12.50 


10.20 


9.10 


7.62 


Bayes Shrink 


37.31 


35.93 


22.40 


15.78 


11.30 


9.22 


8.33 


7.10 


R-Bayes Shrink 


19.64 


19.74 


14.55 


9.35 


7.80 


6.55 


5.75 


5.10 



Other images with it. Table |Vl shows the amount of noise in 
the original image, the image denoised by BayesShrink, and 
the image denoised by R-BayesShrink. To measure the noise, 
we use the variance of a smooth part of the phantom. These 
numbers are the evaluation measure for the amount of noise 
used by radiologists. 

This table shows that R-BayesShrink reduces the noise more 
effectively. Figure [6] shows the CT phantom images denoised 
with R-BayesShrink to compare the edges with the original 
one. All the edges, even the low contrast ones, are kept very 
well. Table |Vll compares the PSNR's of R-BayesShrink and 
BayesShrink. To calculate PSNR, 500mA image is used as the 
reference. This table confirms the strength of R-BayesShrink 
in comparison with BayesShrink. 
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n n □ 



(a) 



(b) 



(c) 



Fig. 6. Phantom images denoised with R-B ayes Shrink. From left to right: 10mA, 50mA, and 300mA (Top: noisy image; Bottom: denoised image by 
R-BayesShrink. 



TABLE VI 

PSNR OF DENOISED CT PHANTOM WITH BAYESSHRINK AND 

R-BayesShrink. 



Tube Current(mA) 


10 


20 


50 


100 


200 


300 


400 


BayesShrink 


14.63 


16.32 


20.18 


22.84 


24.36 


26.52 


26.84 


R-BayesShrink 


16.14 


18.98 


22.31 


24.91 


25.46 


27.60 


27.72 



VI. Conclusion 

We presented a wavelet based denoising method based on 
the structure of the Bayes estimator for GeneraHzed Gaussian 
Distributed (GGD) data. We showed that for this class of 
signals, the Bayes estimate can be well approximated by 
a soft thresholding map. The optimum wavelet coefficient 
soft thresholding map denoted by R-BayesShrink was pro- 
vided. Our proposed threshold value is a generalized form 
of the BayesShrink method. While the threshold value of the 
BayesShrink is only a function of noise variance and noise-free 
data variance, our proposed threshold generalizes BayesShrink 
by considering the shape parameter of the subbands as well. 
Consequently it was shown that BayesShrink is a special case 
of R-BayesShrink when data has a Laplacian distribution. 
Finally we evaluated the proposed method with three different 
data-sets and in all cases the proposed method showed superi- 
ority by simultaneously providing the better PSNR and better 
edge preservation. 
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